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ABSTRACT 

We present the first K'-band, long-basehne interferometric observations of 
the northern Be stars 7 Cas, Per, ( Tau, and k Dra. The measurements were 
made with multiple telescope pairs of the CHARA Array interferometer, and in 
every case the observations indicate that the circumstellar disks of the targets 
are resolved. We fit the interferometric visibilities with predictions from a simple 
disk model that assumes an isothermal gas in Keplerian rotation. We derive fits 
of the four model parameters (disk base density, radial density exponent, disk 
normal inclination, and position angle) for each of the targets. The resulting 
densities are in broad agreement with prior studies of the IR excess flux and the 
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resulting orientations generally agree with those from interferometric Ha and 
continuum polarimetric observations. We find that the angular size of the K' 
disk emission is smaller than that determined for the Ha emission, and we argue 
that the difference is the result of a larger Ha opacity and the relatively larger 
neutral hydrogen fraction with increasing disk radius. All the targets are known 
binaries with faint companions, and we find that companions appear to influence 
the interferometric visibihties in the cases of Per and k Dra. We also present 
contemporaneous observations of the Ha, H7, and Br7 emission lines. Synthetic 
model profiles of these lines that are based on the same disk inclination and 
radial density exponent as derived from the CHARA Array observations match 
the observed emission line strength if the disk base density is reduced by «i 1.7 
dex. 

Subject headings: techniques: interferometric — stars: emission-line. Be — 
stars: individual (7 Cas; 4> Per; C, Tau; n Dra) 

1. Introduction 

Be stars arc rapidly rotating B-type stars that manage to eject gas into a circumstcllar 
disk (observed in H emission lines, an IR excess flux, and linear polarization; Porter & 
Rivinius 2003). Be star disks arc ephemeral and vary on timescales from days to decades 
(in some cases disappearing altogether for extended periods). This inherent time variability 
suggests that gas injection into the disk is only partially the result of the fast spin and 
equatorial extension of the Be star (Townsend, Owocki, & Howarth 2004) and that magnetic, 
pulsational, wind driving, and/or other processes are required for mass loss into the disk 
(Owocki 2005). The evolutionary status of Be stars is still a subject of considerable debate 
(McSwain & Gies 2005; Zorec, Fremat, & Cidale 2005), but there is growing evidence that 
many Be stars were spun up through mass transfer from a binary companion that has since 
become a neutron star, white dwarf, or hot subdwarf (Pols et al. 1991; Gies 2000). 

The direct resolution of Be star disks has come through high resolution radio observa- 
tions {ip Per; Dougherty & Taylor 1992) and long baseline interferometry in a narrow band 
around the Ha emission feature (Stce et al. 2005). The flrst Ha studies focused on the bright, 
northern sky Be stars 7 Cas (Thom, Cranes, & Vakili 1986; Mourard ct al. 1989; Stee et al. 
1995; Berio et al. 1999) and C Tau (Vakih et al. 1998). In a seminal paper, Quirrenbach et al. 
(1997) presented results for seven Be stars based upon narrow-band Ha observations made 
with the Mklll interferometer at Mount Wilson Observatory. They were able to resolve the 
disks around all seven stars, and they found diameters in the range 2.6 — 4.5 mas (Gaussian 
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F WHM) . In most cases the visibilities indicated an elongated shape; for example, in the case 
of C Tau, the major axis was resolved(4.53 mas diameter) but the orthogonal minor axis 
remained unresolved, consistent with the suggestion that we are viewing the disk almost 
cdgc-on. The disk orientation position angle on the sky inferred from polarimetry matched 
the interferometric results in every case. These results have been extended and improved 
through new interferometric observations with the Navy Prototype Optical Interferometer 
(NPOI) in an important series of papers by Tycner and collaborators (Tycner et al. 2003, 
2004, 2005, 2006). 

The relative flux contribution from the disk becomes larger in the infrared spectral 
range due to bound-free and free-free emission from the dense and ionized gas in the disk. 
Waters (1986) developed a disk model for the IR emission that he subsequently apphed to 
determine the disk properties of 101 Be stars based upon their far-IR emission observed with 
the IRAS satellite (Cote & Waters 1987; Waters, Cote, & Lamers 1987), and Dougherty et 
al. (1994) further extended this analysis to study the near-IR excess in 144 Be stars. Howells 
et al. (2001) demonstrated that the near-IR flux excess is correlated with the Ha emission 
strength in Be stars as expected for a common origin in the disk. Thus, the near- and 
mid-lR ranges (Rinehart, Houck, & Smith 1999) offer a new opportunity for resolving the 
disks of Be stars through long baseline interferometry. Stee & Bittar (2001) presented a 
wind model for the infrared emission from Be stars and the corresponding interferometric 
visibihties for the models. Their model for 7 Cas, for example, predicts that the disk should 
appear approximately two times larger in the K-hsaid continuum and in the Br7 emission 
line than found by narrow-band Ha interferometry. Stee (2003) made additional predictions 
for 16 Be stars of the expected X-band interferometric visibilities. The first attempt to 
resolve the infrared disk flux was recently made by Chesneau et al. (2005) who made A^- 
band observations with the VLTl/MlDl instrument of the bright southern Be star a Ara, and 
their results indicate that the disk is smaller than that predicted by Stee (2003), perhaps due 
to disk truncation by a faint binary companion. A'-band interferometric observations from 
VLTI/ VINCI of a Ara were reported by Domiciano de Souza et al. (2003) who interpreted 
the visibilities using a rotationally distorted photospheric model. 

Here we present the flrst near-IR interferometric observations of four bright, northern 
sky Be stars that we obtained with the Georgia State University Center for High Angular 
Resolution Astronomy (CHARA) Array at Mount Wilson Observatory (McAlister et al. 
2005; ten Brummelaar et al. 2005). We describe the interferometric and complementary 
spectroscopic observations in §2 and then we present a simple disk model in §3 that we use 
to predict the interferometric visibflities and spectral line proflles. We discuss the speciflc 
results for each of the four targets in §4. Our results are summarized in §5. 
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2. Observations 
2.1. Targets 

We selected four well known Be stars for this initial program of CHARA Array in- 
terferometry: 7 Cas (HR 264 = HD 5394 = HIP 4427), Per (HR 496 = HD 10516 = 
HIP 8068), C Tau (HR 1910 = HD 37202 = HIP 26451), and k Dra (HR 4787= HD 109387 

= HIP 61281). All the targets except k, Dra have prior narrow-band Ha interferometric 
observations (Quirrenbach et al. 1997; Tycner et al. 2003, 2004, 2006). The first three stars 
were targets in speckle interferometric searches for close companions, and no such compan- 
ions were detected (Mason et al. 1997); however, oscillations in the the proper motion of 
7 Cas suggest the presence of a faint companion with a period of some 60 y (Gontcharov, 
Andronova, & Titov 2000). On the other hand, all four of the targets are known spectro- 
scopic binaries with periods in the range 61 - 204 d. The companion in the case of Per is 
a hot subdwarf whose spectral features appear in the UV spectrum (Gies et al. 1998), but 
for the other three stars the nature of the low mass companion is unknown. 

We summarize the adopted parameters for the target stars in Table 1 that hsts the 
spectral classification, parallax (Perryman 1997), radius, mass, effective temperature, and 
projected rotational velocity (Yang et al. 1990; Abt, Levato, & Grosso 2002) of the bright 
primary star. The next three rows give the radius, mass, and effective temperature of the 
secondary. These are only known for Per and representative values are listed for the others 
assuming that they also have hot subdwarf companions (see §3.3). The next four rows list the 
binary period, epoch of the secondary star's maximum radial velocity (equal to the epoch of 
the secondary's crossing of the ascending node), the primary star to center of mass portion 
of the semimajor axis, and the adopted Roche radius of the primary star. The final row 
provides a key to the references from which these various parameters were adopted. 

2.2. K'-hand Interferometry from the CHARA Array 

The CHARA Array observations were made on various dates during the first two years 
of operation (Table 2). The telescope, instrumentation, and data reduction procedures are 
described in detail by ten Brummelaar et al. (2005). The CHARA Array consists of six 
1 m telescopes in a y-configuration with pairwise baselines ranging from 34 to 330 m in 
length. The Be star observations were primarily made with intermediate to long baseline 
configurations using the K' filter described by McAlister et al. (2005). A summary of each 
night's observations is given in Table 2 that lists the target and calibrator star names, the 
telescope pair (with the maximum baseline given in meters in parentheses; ten Brummelaar 
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et al. 2005), the UT date, and the number of sets of measurements. We made the mea- 
surements using the "CHARA Classic" beam combiner, which is a two-beam, pupil-plane 
(or Michelson) combiner (Sturmann et al. 2003). The fringes were recorded with a near-lR 
detector on either 1 or 2 x 2 pixels at a sampling rate of 100 or 150 Hz depending on the 
seeing conditions. The path length scans were adjusted to obtain approximately five samples 
per fringe spacing. Each measurement set consists of some 200 scans of the fringe pattern 
with photometric calibration scans performed before and after. The raw visibilities were es- 
timated from an edited set of scans using the log-normal power spectrum method described 
by ten Brummelaar et al. (2005). 

We transformed the raw visibilities into absolute visibility V by interleaving the target 
observations with measurements of calibrator stars with small angular diameters (Berger et 
al. 2006). The calibrator stars are generally close in the sky and of comparable K' magnitude 
to the targets, and we list in Table 3 their effective temperatures and gravities as derived 
by other investigators (identified in column 7) . We estimated the angular diameters of the 
calibrator stars by comparing their observed and model fiux distributions. The angular 
diameter of the limb darkened disk 9ld (in units of radians) is found by the inverse-square 
law: 

/a (observed) /d /jn2 1/,2 



F,(emitted) lO-o-^^^ ' ^^^^^ " 4^^^ 

where the ratio of the observed and emitted fluxes (reduced by the effects of interstellar 
extinction Ax) depends on the square of the ratio of stellar radius i?^ to distance d. Most 
of the calibrator stars are relatively close by and the interstellar extinction is negligible. 
The one exception is HD 107193, and in this case we adopted an extinction curve law from 
Fitzpatrick (1999) that is a function of the reddening E{B — V) and the ratio of total-to- 
selective extinction R = Ay / E{B — V) (set at a value of 3.1). 

The model fluxes were interpolated from the grid of models from R. L. Kurucz^ based 
upon the values of Tcfj and \ogg in Table 3. These model fluxes are based upon solar abun- 
dance, plane-parallel, local thermodynamic equilibrium (LTE), line blanketed atmospheres 
with a microturbulent velocity of 4 km s~^ that should be adequate for these A- and F- 
type main sequence stars. We compiled data on the observed fluxes from the optical to the 
near-IR. We found Johnson UBV photometry for all the cahbrators, and these magnitudes 
were transformed to fluxes using the cahbration of Cohna, Bohlin, & Castelli (1996). We 
also included Stromgren photometry for several calibrators using the calibration of Gray 
(1998). Spectrophotometry from Kharitonov, Tereshchenko, & Knyazeva (1988) was used 
for HD 6210. AU of the stars are included in the 2MASS All-Sky Catalog of Point Sources 
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(Cutri et al. 2003), and we converted the JHK magnitudes to fluxes using the cahbration 
of Cohen, Wheaton, & Megeath (2003). In a few cases the errors in the ^AfA S'5' magnitudes 
were unacceptably large, and we relied instead on IR fluxes from the COBE DIRBE instru- 
ment (Smith, Price, & Baker 2004). The angular diameters we derived by fitting the spectral 
energy distributions are given in column 6 of Table 3. Most of these calibrators have angu- 
lar diameters that are small enough that their associated errors will introduce only minor 
systematic errors in our derived target angular sizes (van Belle & van Belle 2005). 

It is normal practice in intcrfcrometry to estimate the actual visibility of a calibrator 
by adopting a uniform disk function for the visibility curve as a function of wavelength and 
baseline, but in fact it is no more difficult to apply a visibility curve for a limb darkened 
disk provided the form of the hmb darkening is known (Davis, Tango, & Booth 2000). 
We used the effective temperatures and gravities from Table 3 to estimate the X-band 
limb darkening coefficients from the tables of Claret (2000), and then we calculated the 
limb darkened visibility curves for the calibrators (for the projected baseline at the time of 
observation and the limb darkened angular diameter given in Table 3) by making a weighted 
sum of the predicted visibilities over the wavelength band of the K' filter (to account for 
bandwidth smearing). We then made a simple linear interpolation in time between the 
target and calibrator observations to estimate the ratio of raw to absolute visibility required 
to transform the results to absolute visibility. Note that we also made a sphne fit of the time 
evolution of calibrator visibility, and the difference between the spline and linear interpolated 
visibilities was added quadratically to the error budget for the final calibrated visibility to 
help estimate the errors introduced by the time interpolation scheme. 

The cahbrated visibilities are presented in Table 4 (available in full in the electronic 
version of the paper). The columns in this table give the target name, the heliocentric Julian 
date of the mid-point of the data set, the binary orbital phase determined from the period 
and epoch given in Table 1, the telescope pair used for the observation, the projected baseline 
and position angle of the target as viewed from the telescope pair at the mid-point time, and 
the calibrated visibility. Note that the visibility errors quoted represent the internal errors 
from the power spectrum analysis of the fringes plus a term introduced for the calibration 
process, and these may underestimate the actual visibility errors. We find, for example, 
that there are some closely spaced observations on certain nights with comparable projected 
baselines and position angles that have a scatter that is several times larger than the formal 
errors. There are 11 such subsets of 4 or more measurements from within a specific night 
that have ranges of 5(basehne) < 10 m and ^(position angle) < 10°, and the average ratio 
of the standard deviation of V within a subset to the mean of the quoted internal error in 
y is a factor of 2.8. Thus, the quoted errors are probably lower limit estimates of the true 
error budget. 
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2.3. Optical Spectroscopy 

The disks of Be stars exhibit large temporal variations in their H-Balmer emission 
strengths that presumably reflect large structural changes in their geometry and/or density 
(Porter & Rivinius 2003). Thus, we obtained new spectroscopic observations of the target 
Be stars in order to compare the Balmer emission levels at times contemporaneous with the 
interferometric observations. The sources and details of the spectroscopy are hsted in Table 5 
that reports the observatory and telescope of origin, the spectroscopic instrument, spectral 
resolving power R = spectral range recorded, UT dates, and the names of targets. 

Most of the observations of the Ha and H7 lines were made with the Kitt Peak National 
Observatory (KPNO) coude feed telescope, and these have moderate resolution and high 
S/N properties. The Ha spectra of k, Dra were obtained at both the University of Toledo 
Ritter Observatory and the Ondfejov Observatory at the same time as the first CHARA 
Array observations were made, and a low dispersion blue spectrum was also obtained at 
that time with the University of Texas McDonald Observatory 2.7 m telescope (recording 
H/3 through the higher Balmer sequence). All of these spectra were reduced using standard 
routines in IRAF^ to create continuum rectified versions of each spectrum. The atmospheric 
telluric lines in the vicinity of Ha were removed from the KPNO coude feed spectra (Gies et 
al. 2002), but these features were quite weak in the k, Dra Ha spectra and were left in place. 
We describe these Balmer line profiles in §4 and §5 below. 



2.4. K-hand Spectroscopy 

The IR excess flux from Be star disks will tend to dilute the photospheric and emission 
lines in the near-IR spectral range, and the strength of the Br7 emission line in particular 
offers a sensitive test of the IR excess fluxes derived from interferometry. Thus, we also 
obtained low resolution A'-band spectroscopy of all the targets except k Dra with the Uni- 
versity of Texas McDonald Observatory 2.7 m telescope and CoolSpec spectrometer (Lester 
et al. 2000). We used a I'.'O entrance slit and a grating with 75 grooves mm~^ in second order, 
yielding a spectral resolving power oi R ^ 1520 (Table 5). The detector was a NICM0S3 
HgCdTe 256 x 256 pixel array. The spectra were reduced using the methods described by 
Likkel et al. (2006) except that in our case the sky background was estimated from the 
intensity at off-star positions along the projected slit. 



^IRAF is distributed by the National Optical Astronomical Observatory, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc. (AURA), under cooperative agreement with the 
National Science Foundation. 
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We also made a number of spectra at different air masses of the broad-lined, AO V stars 
HR 567 (HD 11946) and HR 2398 (HD 46553), which we used to remove the numerous atmo- 
spheric features in this spectral band. The stellar spectra of these two stars are essentially 
featureless with the exception of the photospheric Bv-f line (Wallace & Hinkle 1997), and 
we formed pure atmospheric spectra by dividing our observed spectra by the average stellar 
spectra of the similar stars HR 5793 and HR 7001 (after broadening for the lower resolution 
of our spectra and for the higher projected rotational velocities of HR 567 and HR 2398) 
that we obtained from the spectral library of Wallace & Hinkle (1997)'''. The target spectra 
were then divided by these atmospheric spectra (after small line depth corrections to account 
for differing air masses) to obtain the final continuum rectified versions that are presented 
below (§5). 



The observed interferometric visibilities can be interpreted by assuming that the Be star 
disks have a relatively simple spatial intensity in the sky, perhaps a uniform ellipsoid, ring, 
or Gaussian eUipsoid (Tycner et al. 2006). However, the quality of the interferometric data is 
sufficient to begin to explore more physically motivated models for the disk appearance. Here 
we present a relatively simple model for the near-IR emission that is derived from the Be 
circumstellar disk model of Hummel & Vranckcn (2000). The advantages of this approach 
are that we can compare the resulting disk gas densities with prior work on the near-IR 
spectrum and that we can also predict the disk emission line strengths for comparison with 
our observed spectra. In the following subsections we outline the basic assumptions and 
characteristics of the model that we then apply to fit the CHARA Array observations of 
visibility. We discuss the model predictions for spectroscopy below in §4. 



Hummel & Vrancken (2000) present a simple model for the circumstellar disks of Be 
stars that they use to explore the emission hue shapes. The basic concept is that the disk is 
axisymmetric and begins at the stellar surface. The gas density in the disk is given by 



3. 



Be Star Disk Models 



3.1. Disk Geometry and Model Pctrameters 





''ftp: / / ftp.noao.edu/ catalogs/medresIR/K_band/ 
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where R and Z are the radial and vertical cylindrical coordinates (in units of stellar radii), 
Pq is the base density at the stellar equator, n is a radial density exponent, and H{R) is the 
disk vertical scale height. This gas scale height is given by 

H{R) = ^rI (3) 

where Cs is the sound speed and Vk is the Keplerian velocity at the stellar equator. We 
assume that the outer boundary of the disk occurs at a radius Rd that is equal to the Roche 
radius of the Be star for these four binary targets (Table 1). This assumption is probably 
more critical for the mid-IR emission spectrum, but it not important for the K' emission 
that is generally confined to regions well within the Roche radii of our targets. 

Hummel & Vrancken (2000) assume that the temperature is constant throughout the 
disk and is given by Td — ^T^s where T^s is the stellar effective temperature. In fact, recent 
physical models for Be disks by Carciofi & Bjorkman (2006) show that the gas temperature 
is actually a function of both R and Z, but their work shows that the temperature is about 
60% Tcfr in the outer, optically thin parts of the disk, so we adopt this value for the isothermal 
approximation used here. 

We assume that the star itself is spherical and of uniform brightness in the K' band. 
Neither of these assumptions is probably correct for a rotationally distorted and gravity 
darkened Be star (Townsend et al. 2004), but the photospheres of the targets are small 
enough (see Table 6 below) that these details are inconsequential for the baselines and 
wavelength of the CHARA Array observations. We suppose that the disk is observed with 
the disk normal oriented to our line of sight with an inclination angle i and a position angle 
a measured east from north in the sky. Note that this position angle convention is 90° 
different from that adopted in the Ha interferometric studies where the position angle of the 
projected major axis is usually given. Thus, the four parameters that define a disk model 
are the base density po, the density exponent n, and the orientation angles i and a. 



3.2. K'-band Continuum Images and Interferometric Visibility 

We determine the surface brightness of the disk plus star over a projected rectilinear 
coordinate grid on the sky by solving the equation of transfer along a ray through the center 
of each grid position, 

I^Sd{l-e--) + he-- (4) 

where / is the derived specific intensity, Sd is the source function for the disk gas (taken as 
the Planck function for the disk temperature Td), h is the specific intensity for a uniform 
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disk star (taken as the Planck function for Teff), and r is the integrated optical depth along 
the ray. The main shortcoming in this expression is the neglect of a scattered light term 
to account for Thomson scattering of photospheric intensity (see, for example, cq. [6] in 
Bjorkman & Bjorkman 1994). For the typical electron densities that we find for our Be 
star targets, this term will amount to only a few percent of the stellar specific intensity at 
locations close to the star, which is usually much less than the disk source function in the 
optical thick portions, so this simplification is acceptable in the K' band. 

The disk optical depth in the ncar-lR is predominantly due to bound-free and free-free 
processes, and we can express the increment in optical depth with an incremental step ds 
along a given ray as 

dr = C{X,Td)p{R,Zfds (5) 

where the coefficient C{X,Td) is given by equation (5) in Dougherty et al. (1994) (see also 
Waters 1986). This coefficient includes terms for the Gaunt factors for bound- free and free- 
free emission that we evaluated for the K' band using the tables in Waters & Lamers (1984). 
We adopted an ionization model for the disks of the two hotter Be stars, 7 Gas and Per, 
that assumed ionized H, singly-ionized He, and doubly-ionized G, N, and O atoms, while 
for the cooler Be stars, ( Tau and k, Dra, wc assumed ionized H, neutral He, and singly 
ionized C, N, and O. Note that we evaluated the optical depth coefficient only at the central 
wavelength of the K' filter since the coefficient varies slowly with wavelength. Also note 
that we have accounted only for continuum emission in this band since the Br7 emission 
contribution is small compared to the fiux integrated over the K' band (§5). 

The image of the Be star in the sky in constructed by integrating the optical depth along 
each ray according to equations (2) and (5), and then populating each pixel of the image 
using equation (4). The pixel scale on the sky is set by the adopted angular diameter of the 
star (derived from the parallax and radius in Table 1). Wc show an example calculation for 
7 Gas in the left panel of Figure 1. Portions of the disk projected against the sky include only 
the first term of equation (4), and the disk intensity attains a maximum of Sd ~ 55%/* in the 
inner optically thick regions. At a radius where the optical depth is approximately unity, the 
disk becomes optically thin and fades with increasing radius. Part of the background disk is 
occulted by the star while the part in foreground attenuates the photospheric intensity. We 
also include a binary star option to add the disk of a secondary into the image at a position 
appropriate for the time of observation. 

The interferometric visibilities are then calculated by making Fourier transforms of such 
disk images, and we did this following the techniques described by Aufdenberg et al. (2006). 
The comparison with the observed visibility is done by determining the discrete Fourier 
transform of the image for the {u, v) spatial frequency pair of the observation (see eq. [19] 
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in Aufdenberg et al. 2006), and this calculation is done successively and summed over five 
wavelength bins according to the transmission curve of the K' filter (McAlister et al. 2005). 
This procedure assumes that spectral energy distributions in the K' band of our targets are 
similar to the stellar transmission curve (Rayleigh- Jeans tail), and this is probably acceptable 
in the K' band where the disk emission is just beginning to become an important component 
of the total flux (Waters 1986). The Fourier transforms along the major and minor axes of 
the projected disk image for the 7 Cas example are illustrated in the right hand panel of 
Figure 1 (here without accounting for bandwidth smearing) . The top curve shows the slower 
decline in V for the smaller minor axis while the bottom curve illustrates the rapid decline 
for the better resolved major axis. Note that the star itself is so small {Olb — 0.50 mas) 
that its visibility curve in the absence of a disk would decline to only V — 0.87 at a baseline 
of 300 m. 

Figure 2 shows how the image and visibilities change by altering the disk inclination 
angle from 51° to 80° (the dotted lines in the right hand panel copy the original visibility 
curves from Fig. 1). The rays through the outer positions along the major axis now traverse 
more optical depth because of the oblique view through the disk, and this brightens the outer 
regions making the disk effectively larger in this dimension (so that the decline in V with 
baseline is steeper). On the other hand, the minor axis now appears more foreshortened and 
the visibility decline is reduced. Thus, the V variation with position angle in the sky (with 
respect to the major axis) is very sensitive to the disk inclination angle i. 

The changes resulting from a decrease in disk base density po ^"^^ illustrated in Figure 3. 
Now the density is lower throughout the disk, and consequently the optically thick part of 
the disk appears smaller in both dimensions (and the V curves decline more slowly with 
baseline). Changes in the opposite sense result from a decrease in the density exponent as 
shown in Figure 4 (from n — 2.7 to 2.0). The radial density decline is less steep in this case 
so that much more of the inner disk is optically thick, leading to a larger appearance in the 
sky (steeper decline in the V curves). These examples show that the visibility curves are 
most sensitive to the disk size as given by the boundary between the optically thick and thin 
regions. Since the optical depth unity radius in the disk plane depends on pQr~^"+'^/^, there 
will be a locus of parameters that will produce visibility curves with similar decline rates in 
the main lobe and with subtle differences appearing only at intermediate and longer baselines. 
Thus, accurate determinations of both po and n will require visibility measurements over a 
broad range of baselines. 
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3.3. Model Fitting Results 

We found best fit solutions for the four disk parameters by fitting the model and observed 
K' visibilities using a numerical grid search technique. We used an iterative search on each 
of the individual parameters to find the set that provided the lowest overall statistic. 
The results are fisted in Table 6 for both single star and binary star models (see below). 
The seventh row of Table 6 reports the minimum values, and these are all much larger 
than the expected value of unity for acceptable fits. Recall that the internal errors on V may 
underestimate the actual errors by a factor of 2.8 for these measurements (§2.2), so we might 
instead expect to find best fits with xl ~ 2.8^ = 7.84. This is more or less consistent with 
the fitting results for all but 7 Cas, where the residuals are still larger than expected. We 
estimated the errors in the fitting parameters by renormalizing the minimum x^ to unity and 
then finding the excursion in the parameter that increased to a value of 1 + 1/ [Ny — 4) 
(for Nv visibility measurements and 4 fitting parameters) . Since the disk parameters po ^-nd 
n are so closely coupled in the fits (§3.2), we searched the locus of the pair- wise fits for the 
same increase in x^, and we quote the half-range in these parameters for their errors. 

The first two rows of Table 6 give the orientation parameters, the disk plane normal 
inclination i and position angle a (east from north). Note that there is a 180° ambiguity 
in the value of a (i.e., whether the disk north pole is in the eastern or western half of the 
sky), so we assume the range < a < 180°. The third and fifth rows list the logarithm 
(base 10) of the disk base density log po and the disk density exponent n, respectively (base 
electron density is logA^e = logPo + 23.6). These are compared with the same parameters in 
rows four and six, respectively, estimated from model fits of the IR fiux excess from IRAS 
observations that were computed by Waters et al. (1987). Waters et al. (1987) make a 
number of simplifying assumptions about the disk geometry that make a direct comparison 
of results unreliable, but in general there is reasonable agreement except in the two cases 
where a binary companion may influence the results, Per and k Dra (see below). Row 8 
gives a derived model quantity related to the IR flux excess, E(y — K) — 2.5 log(l + F^/F^), 
where F^ and F^, are the integrated disk and stellar flux, respectively. If we assume that 
the disk fiux is relatively small in the optical, then this quantity can be directly compared 
to the observed IR fiux excess measurements by Dougherty et al. (1994) that are listed in 
row 9. The model and observed values of E{V — K) are comparable but not always a good 
match. We suspect the differences may relate to the temporal variations in the disk emission 
(from ~ 1990 for the IR magnitudes to ~ 2005 for the CHARA Array measurements). Rows 
10 and 11 compare the adopted stellar angular diameter Old with the model disk emission 
FWHM diameter along the projected major axis. The disk emission FWHM was estimated 
from the width at 50% of the emission intensity in the optically thick part just above the 
photosphere. 
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We know that all four targets are close binaries with faint, low-mass companions, but 
the nature of the companion is only known in the case of Per (Gies et al. 1998). The 
secondary of Per is a hot helium star, the stripped down remains of what was once the 
more massive star before mass transfer began. If we assume that all our Be star targets arc 
likewise the mass gainers that have been spun up by the angular momentum carried with 
mass transfer, then they may also have helium star companions (or possibly white dwarf 
or neutron star companions). We present here a sample calculation of how such a binary 
companion would affect the K' visibilities. For the purpose of this demonstration, we assume 
that the companion is a hot (Teff = 30000 K) and small {R — IRq) helium star (except in the 
case of (f) Per where we have better estimates for these parameters; see Table 1). According 
to the interacting binary scenario, the angular momentum vector of the mass gainer Be star 
should be parallel with the orbital angular momentum vector, so the disk plane corresponds 
to the orbital plane. Consequently, we can apply the orientation parameters for the disk, i 
and a, also to the orbit of the secondary. We know the epoch of the secondary's ascending 
node passage from the spectroscopic radial velocity curves (Table 1), so the only remaining 
unknown for the secondary's astrometric orbit is the sense of motion on the sky (increasing 
or decreasing position angle with time). We follow the normal convention of assigning an 
inclination in the range 0° < i < 90° for counterclockwise apparent motion and in the range 
90° < i < 180° for clockwise motion (so that the longitude of the node is Q = a + 270° for 
i < 90° and fl = a + 90° for i > 90°). We made fits of the visibility for these two binary orbit 
cases, and we show in Table 6 the results for the best fit sense of motion. The minimum xl 
was significantly reduced only for the case of Per (where the properties of the secondary 
are known; the estimated magnitude difference is AK' = 2.9 mag). The inclusion of the 
binary companion also resulted in a significant change in the fitted disk properties for k Dra 
(discussed further in §5.4). 

It is valuable to compare the K' angular diameters with those derived from narrow- 
band Ho; interferometry, and to do so we need to fit the visibilities with a simple Gaussian 
ellipsoidal model as done for the Ha results. We made these fits also by a numerical grid 
search of the four parameters that describe the Gaussian eUipsoidal model: the angular 
FWHM of the major axis ^gd, the ratio of the minor to major axes r (> cosi), the position 
angle of the major axis (j) {— a ± 90°), and the fraction of the total fiux contributed by 
the photosphere of the Be star cp. The best fit results are listed in Table 7 for these 
single star models. Our results are compared with those from Quirrenbach et al. (1997) and 
Tycner et al. (2004, 2006) based on interferometric observations of the Ha emission, and the 
orientation parameters appear to be in reasonable agreement. We find, however, that the 
K' disk diameters are somewhat smaller than the Ha diameters 64% as large), contrary 
to the wind model predictions of Stee & Bittar (2001). Note that the xl residuals of the 
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Gaussian elliptical fits are comparable to those from the physical model fits (Table 6), and 
additional longer baseline observations will be required to distinguish between them. 



We can also use the simple disk model to estimate the line profiles of the disk emission 
lines (Hummel & Vrancken 2000). However, the problem is now more complex because of the 
wavelength dependence of the emission, and a number of key assumptions about the atomic 
populations must be made. Despite these difficulties, it is nevertheless interesting to check 
whether or not the disk density parameters derived from interferometry lead to predicted 
line profiles that are consistent with the observed ones. Here we discuss the line synthesis 
predictions for the disk models and their comparison with our spectroscopic observations. 



The line synthesis procedure has many similarities with the disk image calculations 
present above. We assume the same basic density law (eq. [2]), the same creation of a 
rectilinear grid of points covering the projected positions of the star and disk in the sky, 
and the solution of the equation of transfer (eq. [4]) along along a ray through the center of 
each grid point. However, for the line synthesis the expression for the optical depth (eq. [5]) 
becomes (Hummel & Vrancken 2000) 



and now we consider a wavelength dependent optical depth t\ for the region surrounding 
an emission line (we adopt an equivalent radial velocity grid in 10 km s~^ increments from 
—2000 to +2000 km s^^). The first term in the expression is the scattering coefficient for a 
classical oscillator that is multiplied by the oscillator strength / and the rest wavelength Aq 
for the specific transition (H7, Ha, and Br7). 

The next term N[Hi) is the local number density of neutral H atoms in energy level i. 
This can be set approximately by balancing the recombination rate with the photoionization 
rate (for the case when all the bound-bound transitions are very optically thick; Cassinelli, 
Nordsieck, & Murison 1987) 



4. Model Emission Lines and Spectroscopy 



4.1. 



Disk Emission Line Models 



dTx = /Ao N{Hi) Px{Vr,N,) ds 



(6) 



00 




(7) 
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where ck{T) is the recombination coefficient, Ne and Np are the electron and proton number 
densities, Ni is the number density of hydrogen in energy level i, is the photoionization 
cross section, and Jy is the mean radiation field. If we assume that the latter results only 
from unattenuated starlight, then the mean field at a distance r from the star can be written 

J. = W{t)II (8) 

where W{r) is the dilution factor given by 

iy(r) - 0.5(1 - [1 -r-2]5), (9) 

r is the radial distance from the star r = {E? + ^^)^, and J* is the stellar specific intensity. 
Cassinelli et al. (1987) rewrite the photoionization/recombination equilibrium to express the 
level population Ni in terms of a variable g^, 

oo 

TV. = N,Np a{T)/ j ^a,W{r)It dv = N^N^q,. (10) 

In the stellar photosphere, l^(r) = 1, /* = Bi,{Tes) (the Planck function), and the gas will 
be in LTE, so that q^^^ — N^/ (NgNp) will be given by the Boltzmann-Saha equation. In the 
circumstcllar disk, the photoionization rate will vary with W{r) and the recombination rate 
will depend on a{T) oc T~^/^, so level populations will be given by 



N^ = g[™V^efr/T^ N^Np/Wir). (11) 

The actual populations probably have a more complex dependence on position in the disk 
(Carciofi & Bjorkman 2006). 

The final term in optical depth equation is the normalized line profile Px{Vr, Ne) that 
depends on the local radial velocity Vr, electron density, and the adopted temperature of 
the isothermal disk T^. We adopted the convolved thermal and Stark broadened profiles 
calculated by T. Schoning and K. Butler^ for Ho; and H7, but we used only a simple ther- 
mal profile for Br7. These profiles were collected into a matrix of 401 wavelength points 
(corresponding to our radial velocity grid), 60 positive Doppler shifts (in increments of 10 
km s~^), and four electron number densities (logiVe = [10.5, 12, 13, 14]). 

The line synthesis calculation proceeds by solving for the intensity spectrum from each 
pixel in the projected image and then summing these to create a predicted fiux profile. 



^http: / / nova.astro.unid.edu/Synspec43/data/hydprf.dat 
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At each position along a ray through a pixel, we determine the Keplerian velocity and 
Doppler shift, gas density, H population, and profile shape (simply adopting the closest 
match in the grid), and then we integrate the optical depth for each wavelength point. The 
summed spectral intensity is then found through equation (4), but in this calculation J+ 
is replaced with a specific intensity corresponding to a photospheric absorption line with 
the corresponding value of surface normal orientation and Doppler shift for the projected 
position. These intensity profiles were calculated by first creating fine blanketed LTE model 
atmospheres using the ATLAS9 program written by Robert Kurucz. These models assume 
a stellar temperature and gravity derived from the parameters in Table 1, solar abundances, 
and a microturbulent velocity of 2 km s~^. Then a series of specific intensity spectra were 
calculated using the code SYNSPEC (Hubeny & Lanz 1995) for each of these models and 
H lines. These profiles include the effects of limb darkening for the stellar photosphere. 
Because our model assumes axial symmetry, we complete the calculation for only the half 
of the projected image with positive radial velocity and then we add a mirror-symmetric 
version of the fiux profile to account for the other half with negative radial velocity. The 
resulting profiles are then rectified to unit continuum fiux and smoothed by convolution with 
a Gaussian function to match the instrumental broadening of our spectroscopic observations 
(Table 5). We also need to account for the disk IR continuum excess in the case of Br7, so 
we renormalize the continuum fiux Fx to 

F', = {F, + e)/{l + e) (12) 
where e = iQO-mv-K) _ ^ (^(y _ ^) gj^g^^ rp^^^g 

4.2. Emission Line Fitting Results 

The model emission profiles depend on three parameters, the disk base density po, the 
density exponent n, and the disk inclination i. Since our main interest is in checking the 
consistency of the interferometric and spectroscopic results, we set the values of n and i from 
the fits of the CHARA Array visibilities and then varied po to find the best match with the 
observed line profiles. Our results are gathered in Table 8 for each of the four targets. The 
first three rows give the measured equivalent widths (where a negative sign indicates a net line 
flux above the continuum) for Ha, H7, and Br7, respectively. The next three rows give the 
corresponding disk base density po that yields synthetic proflles with the observed equivalent 
width (or that best match the profile structure in the case of H7 where the equivalent width 
is complicated by the presence of line blends) . The final row lists for comparison the value of 
Po as derived from the interferometric data (single star results for 7 Gas and C Tau; binary 
star results for (f) Per and k, Dra). We see that there is relatively good agreement between 
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the densities derived from the three different H hnes, and specifically the near agreement of 
the result from Br7 with those from the Balmer lines supports the flux renormalization near 
Br7 (eq. [12]) set by the interferometric determination of the IR excess. The emission line 
and interferometric densities vary in tandem among the analyses of the four stars, but we 
find that the emission line densities are systematically lower by approximately 1.7 dex. We 
suspect that the problem lies in the simplying assumptions made in the model concerning 
the isothermal disk temperature and H populations that both probably lead to overestimates 
of the line emission from the outer parts of the disk. 

The synthetic model profiles are directly compared with the observed ones in the next 
section, and we find that all the strong emission lines (formed in optically thick regions) 
appear to be much broader than predicted. We show one example of this difference in 
Figure 5 for the Ha line in 7 Cas. The basic synthetic model that matches the observed 
equivalent width shows the strong double peaks associated with Keplerian motion in the 
disk. The observed profile, on the other hand, lacks well defined peaks and shows much 
broader line wings. We also show in Figure 5 a smoothed version of the synthetic profile 
formed by convolution with a Voigt function fit using a damping parameter F = 934 km s~^, 
and this broadened version matches the observed profile quite well. Thus, it appears that 
the line broadening due to Keplerian motion in the disk is insufficient to explain the shapes 
of the strong lines (although the weaker H7 profiles are in better agreement; see §5). There 
are a number of possible explanations for the discrepancy: 

(1) We have assumed that the Ha source function is constant with radius in the disk, but more 
detailed models by Carciofi & Bjorkman (2006) show that the temperature and presumably 
the source function decrease with radius. This suggests that the Ha emission is probably 
stronger in the inner disk than portrayed in our model, and since large orbital motions in 
the inner disk produce the line wings, the model almost certainly underestimates the wing 
emission. 

(2) Macroscopic motions such as stellar wind outfiow (Stee et al. 1995) and turbulence may 
be important (especially close to the star where the emission is particularly strong). 

(3) Non-coherent scattering of photons into the wings of strong emission lines may extend 
their widths by as much as a few hundred km s~^ (Hummel & Dachs 1992; Hummel & 
Vrancken 1995), but this may be insufficient to explain the extreme broadening observed in 
Ha. 

(4) Thomson scattering by free electrons in the inner disk may create extended line wings 
(Poeckert & Marlborough 1979). 

(5) CoUisional Stark broadening will produce similar extended wings, but tests with our 
model suggest that the base density would need to be increased by about two orders of 
magnitude for the Stark wings to become prominent. This might be possible if disk gas 
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consists of dense clumps. 

(6) Raman scattering of photons from the vicinity of the ultraviolet Ly/3 line could lead to 
broad wings in Ha (Yoo, Bak, & Lee 2002) if there is a sufficient column density of neutral 
hydrogen in the disk (A^(H 1) ~ ]^q20-23 Q^a'^). This process might be especially applicable 
to edge-on and/or cooler disks. 

5. Individual Tcirgets 
5.1. 7 Cas 

The star 7 Cas is perhaps the best documented Be star in the sky, and it has displayed 
Be, Be-shell, and normal B-type spectra over the course of the last century (Doazan 1982). 
It currently displays strong Ha emission, and the Ha disk size was recently determined 
through observations from NPOI (Tycner et al. 2006). Our best fit disk (Table 6) and 
Gaussian elliptical models (Table 7) are generally consistent with the Ha interferometric 
results and indicate that we view this Be star disk at an intermediate inclination. The disk 
density exponent we derive, n = 2.7 ± 0.3, matches within errors the value n = 2.8 ± 0.1 
found by Hony et al. (2000) based on an analysis of the 2.4 - 45 fim spectrum from ISO- 
SWS. We show the model K' image of 7 Cas in the left hand panel of Figure 6 and the 
corresponding visibility image (without the effects of bandwidth smearing) is shown in the 
right hand panel. Both images are scaled logarithmic representations made to accentuate the 
low level features. The {u, v) frequencies of the CHARA Array observations are shown in the 
upper half of the visibility image with symbol sizes proportional to the observed visibility. 
They also appear in point-symmetric locations in the lower half where they represent the 
residual-to-error ratio (vertical segment for Vobs > Vcai and horizontal for Vohs < Kai)- The 
residuals from the fit are large for our observations of 7 Cas, but they do not appear to 
indicate any systematic deviation from the model since we find some of the largest positive 
and negative deviations in the same part of the (u, v) diagram. We suspect that the high 
value of xl refiects an incomplete accounting of the full error budget. 

Harmanec et al. (2000) discovered that 7 Cas is a single-lined spectroscopic binary, and 
their orbit was confirmed in subsequent Doppler shift measurements by Miroshnichenko et al. 
(2002). We made a fit of the visibilities by including a possible hot subdwarf companion that 
was placed in the model images using the spectroscopic ephemeris from Miroshnichenko et 
al. (2002) (Table 1). This modification made a shght improvement in the residuals of the fit, 
and models with a secondary hotter than T^s — 30 kK showed somewhat better agreement. 

The profiles of the Ha, H7, and Br7 emission lines are plotted as thick lines in Figure 7. 
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We also plot the synthetic model profiles as thin lines for each feature. The model profiles 

of the weaker H7 and Br7 features make a reasonable match of the observed profiles, but, 
as noted above, the Ha profile is much broader and smoother than predicted. Robinson, 
Smith, & Henry (2002) suggest that there may be significant turbulence in the inner disk of 
7 Cas that is associated with magnetic dynamo processes, and such turbulent motions could 
help to explain the wings observed in Ha. 

5.2. Per 

The binary system Per is the only one of the four targets where we known the nature of 
the secondary and can reliably predict its flux contribution in the near-lR (Gics ct al. 1998). 
The companion is a hot subdwarf whose spectrum appears clearly in UV spectroscopy from 
HST, and from the relative strengths of its spectral lines, we can estimate with confidence its 
temperature, radius, and hence its near-IR flux contribution. Indeed the fit of the CHARA 
Array visibihties is better for the binary model than the single star model (Table 6). Here 
the inclusion of the binary flux leads to a downwards revision of the disk size. We show 
in left hand panel of Figure 8 the image of the entire binary for the observation made on 
HJD 2,453,656.9 (orbital phase 0.644), which is similar to the binary orientation for most 
of the observations (Table 4). The right hand panel shows the visibility image that is now 
modulated with the interference pattern introduced by the binary companion. Since the 
presence of the secondary does affect the visibility pattern, additional interferometric data 
should lead to an astrometric orbit that will allow us to test whether or not the disk and 
orbit are co- planar (Clarke & Bjorkman 1998). 

The position angle of the disk that we find is about 20° different than that from the 
recent narrow-band Ha measurements from Tycner ct al. (2006) (Table 7). Our solution is 
somewhat influenced by our assumptions about the properties of the secondary's orbit, and 
additional observations at other binary phases will be required to settle the issue of the disk 
and orbit orientations. We also find a much lower density exponent n — than the value 
oln — 3.0 found by Waters et al. (1987) based on fitting the IR-excess from IRAS data, but 
their result is sensitive to the adopted disk radius. For example. Waters (1986) demonstrated 
that a density exponent oin = 2.4 also fit the IR-excess provided the disk radius was limited 
to = 6.5i?0, so the difference in our results probably hinges on differences in the model 
assumptions. 

Figure 9 shows the Ha, H7, and Br7 profiles for (f) Per. The peaks in the H7 profile 
reversed in strength between our observations in 2004 and 2005. This is an important 
reminder that the disk is time variable and is also probably somewhat asymmetrical in 
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appearance (Hummel & Vrancken 1995; Hummel & Stefl 2001), properties that are beyond 
the scope of our simple disk models. Both the Ha and Br7 line are observed to be broader 
than predicted by the disk model. 

5.3. C Tau 

Harmanec (1984) derived the best radial velocity curve to date for this single-hned 
binary, and we adopt his period here but we took the epoch from the the more recent 
spectroscopic study by Kaye & Gies (1997). A comparison of the single star and binary 
fits to the interferometry shows essentially no differences in the resulting xl values, so we 
will discuss the single star results here. The spatial and visibility images resulting from the 
fit are shown in Figure 10. Our results indicate that wc arc viewing this system close to 
edge-on, and the decline in visibility is especially marked along the direction of the projected 
major axis of the disk. Our derived orientation in the sky is in reasonable agreement with 
the Ha results from Tycner et al. (2004), and the density law agrees within errors with that 
obtained from fitting the IR-excess fiux (Waters et al. 1987). 

The spectral line profiles of ( Tau are illustrated in Figure 11. Again we find that 
the observed Ha line is much broader than the model profile, and here we see evidence of 
a significant asymmetry in the disk in the differing heights of the blue and red emission 
peaks. The long term variations in these asymmetries are documented by Guo et al. (1995). 
Vakili et al. (1998) made Ha interferometric observations in narrow bands matching the 
blue and red peaks, and their results support the idea that the disk density is modulated by 
a prograde one-armed oscillation. The H7 line shows a deep central "shell" feature that is 
consistent with the formation of such shell lines in Be stars observed near i = 90° (Hanuschik 
1996). The fact that our model H7 line depth is shallower than the observed depth probably 
indicates that the outer parts of the disk are cooler than assumed in our isothermal model. 

5.4. K Dra 

Juza et al. (1991) and Saad et al. (2005) show that k Dra is a low amplitude, single- 
lined spectroscopic binary with a period of 65.6 d. Saad et al. (2004) describe a long-term 
cyclic variation of about 22 y that is present in the emission lines, and they suggest that this 
period may correspond to the disk's modulation due to an interaction with a third star in an 
eccentric orbit. The semimajor axis of the third star would be on the order of 100 mas for 
the adopted distance, but we find no obvious evidence (separated fringe packets) for such a 
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companion even at the shortest basehnes. Thus, we have confined our test fits here to the 
single star and binary star (P = 65 d) cases. 

The fit for the single star case is odd in several respects. The density exponent is 
unusually small and the inclination oi i — 90° is at odds with estimates from spectroscopy 
{i — 35° — 45°; Juza et al. 1991). The binary fit has similar residuals, but it brings the 
density exponent and inclination into better accord with expectations. Thus, wc will adopt 
the binary model that is based on a hot, subdwarf companion, but we caution that additional 
interferomctric observations are required to make reliable estimates for the companion's orbit 
and flux contribution. 

The model spatial and visibility images arc shown in Figure 12, and we see in this case 
that the binary modulation is equal to or greater than the disk variation in importance. The 
CHARA Array measurements presented here are the first interferomctric observations of this 
target, so there are no Ha results for comparison. Quirrenbach et al. (1997) found that the 
intrinsic continuum polarization angle is generally perpendicular to the projected disk major 
axis in the sky. Clarke (1990) estimated the intrinsic polarization angle for k, Dra as 23° east 
from north (in agreement with the estimate of 27° from a recent polarization analysis by K. 
Bjorkman, private communication), and this agrees within errors with the position angle of 
a = 21° ± 3° from the fit of the interferomctric visibilities. 

The line profiles are shown in Figure 13. We find that the Ha emission strength is 
approximately that expected for the epoch of our observations and the long term variation 
(see Fig. 4 of Saad et al. 2004). Once again the Ha profile (and perhaps the H7 profile) 
appear to be broader than predicted by the Keplerian model. 

6. Conclusions 

The CHARA Array observations represent the first K' interferomctric measurements of 
northern Be stars, and we find that the disks of all four targets are at least partially resolved. 
The derived disk inclination angle and position angle in the sky generally agree with similar 
estimates from narrow-band Ha interferometry where available. We find, for example, that 
the disk of ( Tau appears with a nearly edge-on orientation, which is consistent with the 
prediction for Be stars with shell spectral features (Hanuschik 1996). All four of the targets 
are close binaries with faint companions, and these companions appear to influence the K' 
visibilities of (f) Per and k Dra, an interesting but complicating factor for the interpretation 
of the disk properties. We find, for example, that the radial gradient in disk density (density 
exponent n in Table 6) is smaller in binaries with smaller semimajor axes, which suggests 
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that binary companions do influence the disk properties. 

We developed a simple Keplerian disk model to create synthetic angular images of 
the disks, and we used these models to create model interferometric visibility patterns and 
spectral emission line profiles for comparison with our observations. These models assume 
isothermal disks and ignore scattering in the radiative transfer process, two significant lim- 
itations that should be addressed in a more complete model. The disks are characterized 
by two physical parameters, the base density and density exponent, and two observational 
parameters, the disk normal inclination and position angle. We fit the CHARA Array visi- 
bilities for all four parameters for both single and binary star models, and our results are in 
good agreement with prior studies of the IR-excess in these stars. The parameters derived 
from the interferometry lead to model emission lines that match the observed ones if the 
base densities are reduced by a factor of about 1.7 dex. This mismatch is not unexpected 
considering the many simplifying assumptions made in the model. 

We find that the angular diameters of the Be disks in the X'-band are consistently 
smaller than those found from Ha interferometry. We think that this difference results from 
the larger Ha opacity and the hydrogen ionization structure in the disk. The K' disk emission 
comes primarily from free-free processes in the mainly ionized gas, and the free-free optical 
depth will vary with disk radius according to the square of the density. In the thin disk 
approximation, the optical depth at radial distance R will vary as r oc pp /^-2n-+3/2y' ^^g^ 
(see eq. [2] and [3]). On the other hand, the Ha emission depends on the relatively low 
neutral hydrogen number density, and we assumed above that this varies inversely with the 
stellar dilution factor (eq. [9]; W{r-) oc R^^ at large radius and Z = 0). Thus, we expect the 
Ha optical depth to vary approximately as r oc Po i^^^^+^/^+^Z cosi, i.e., with a shallower 
effective density exponent of n — 1. Consequently, the Ha optical depth reduction with 
increasing disk radius will be less pronounced than in the case of the IR optical depth, 
resulting in a larger spatial extension in the Ha emission compared to that for the K' 
emission. 

This expectation is borne out in the predicted spatial images. We show an example of 
the derived spatial structure in Figure 14 for our single star model of 7 Cas. The sohd fine 
shows the summed K'-hand intensity of the image projected onto the major axis of the disk 
and measured in stellar radii from the center of the star. We used our model to create an Ha 
image adopting the base density from the fit of the observed emission line (Table 8). The 
intensity was formed by summing across a 2.8 nm wavelength band centered on Ha in order 
to match the Ha interferometric observations of Tycner et al. (2006). Note that changing the 
adopted bandwidth will simply result in a rescaling of the contrast between the central star 
and the surrounding disk emission. The disk image appears darker along the line of nodes 
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and along the line of sight because of the minimal shear broadening that occurs at those 
disk azimuthal angles and the subsequent reduction in the range of wavelength over which 
Ha is optically thick. The dashed line in Figure 14 shows the summed Ha image intensity 
projected onto the major axis. We see that Ha emission is indeed more extended along the 
major axis and that the half maximum intensity radii are located at 2.6R^ and TAR^, for 
the K' band and Ha summed intensities, respectively. The ratio of these half maximum 
intensity radii is 0.35, which is comparable to the observed ratio of 0.55 of the FWHM of 
the major axes of the X'-band and Ha emission in the Gaussian eUiptical fits (Table 7). 
Thus, the differences in the radial exponent dependence of density squared and the neutral 
hydrogen fraction may partially explain the larger appearance of the Ha disks. 

The CHARA Array observations demonstrate the promise of these kinds of measure- 
ments for our understanding of Be stars. Many Be stars were probably spun up by mass 
transfer in a past close interaction, and future observations over a range in {u, v) and orbital 
phase coverage will reveal the nature of their companions. The disks of Be stars are driven by 
dynamic processes and evolve over timescales of months to years, so long-term observations 
will help elucidate the physics of disk formation and dispersal. There is ample evidence that 
the disks are structured and prone to the formation of spiral arm features (Okazaki et al. 
2002). Furthermore, in Be binaries like Per, the portion of the disk closest to the hot 
companion may appear brighter due to local heating by the companion's radiation (Gies et 
al. 1998; Hummel & Stefl 2001). The techniques of phase closure using multiple telescopes 
in long baseline interferometry (Kraus et al. 2005) will soon give us the means to study these 
disk processes directly among the nearby Be stars. 
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been provided from the GSU College of Arts and Sciences and from the Research Program 
Enhancement fund of the Board of Regents of the University System of Georgia, adminis- 
tered through the GSU Office of the Vice President for Research. We gratefully acknowledge 
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Table 1. Adopted Stellar Parameters 



Parameter 


7 Cas 


(p rer 


C lau 


K Dra 


Spectral Classification 


B0.5 Ve 


B0.5 111-Ve 


B2 lllpe 


B5-6 llle 


71 (mas) 


5.32 


4.55 


7.82 


6.55 


Ri (Rq) 


10.0 


7.0 


5.5 


6.4 


Ml (Mo) 


15.5 


9.3 


11.2 


4.8 


Ti eff (K) 


28840 


29300 


19000 


14000 


Vsini (km s~^) 


295 


410 


320 


170 


R2 {RqT 


1.0 


1.33 


1.0 


1.0 


M2 {MqY 


1.0 


1.1 


1.0 


0.8 


T2 eS {Kf 


30000 


53000 


30000 


30000 


^ (d) 


203.59 


126.67 


132.91 


61.56 


Ts+ (HJD-2,400,000) . 


50654.3 


50155.1 


46417.3 


50011.0 


Oi sini {Rq) 


15.279 


24.942 


27.562 


10.078 


Rd {R&) 


214 


130 


146 


67 


References 


1,2 


3 


4, 5,6 


7 



^Example parameters for an assumed subdwarf companion for all but 
(f) Per where the estimates are from Gies et al. (1998). 



References. — 1. Harmanec et al. (2000); 2. Miroshnichenko et al. 
(2002); 3. Gies et al. (1998); 4. Harmanec (1984); 5. Yang et al. (1990); 
6. Kaye & Gies (1997); 7. Saad et al. (2005). 
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Table 2. CHARA Array Observations 



Target 


Galibrator 


Baseline 


Date 




No. of 


Name 


Name 


(m) 


(UT) 




Sets 


7 Gas 


HD 6210 


W1/S2 (249) 


2003 


Oct 


07 


5 


7 Gas 


HD 6210 


W1/S2 (249) 


2003 


Oct 


10 


6 


7 Gas 


HD 6210 


W1/S2 (249) 


2003 


Oct 


14 


5 


7 Gas 


HD 6210 


W1/S2 (249) 


2003 


Oct 


16 


6 


7 Gas 


HD 6210 


W1/W2 (111) 


2005 


Oct 


12 


12 


7 Gas 


HD 6210 


W1/W2 (111) 


2005 


Oct 


13 


11 


7 Gas 


HD 6210 


W1/W2 (111) 


2005 


Oct 


15 


11 


Per 


HD 6961 


W1/S2 (249) 


2003 


Oct 


03 


5 


Per 


HD 6961 


W1/W2 (111) 


2003 


Oct 


12 


10 


(f) Per 


HD 6961 


W1/W2 (111) 


2003 


Oct 


13 


7 


Per 


HD 6961 


W1/W2 (111) 


2003 


Oct 


14 


14 


C Tau 


HD 43042 


Sl/El (331) 


2004 


Dec 


21 


6 


CTau 


HD 32977 


Sl/El (331) 


2004 


Dec 


22 


1 


( Tau 


HD 43042 


Wl/Sl (279) 


2005 


Apr 


01 


2 


C Tau 


HD 43042 


W1/W2 (111) 


2005 


Dec 


04 


8 


C Tau 


HD 43042 


W1/W2 (111) 


2005 


Dec 


05 


9 


K Dra 


HD 107193 


S1/S2 (34) 


2005 


Apr 


02 


13 


K Dra 


HD 107193 


Sl/El (331) 


2005 


Apr 


06 


5 


K Dra 


HD 107193 


Sl/El (331) 


2005 


Apr 


13 


1 


K Dra 


HD 107193 


W1/W2 (111) 


2005 


Dec 


04 


3 


K Dra 


HD 107193 


W1/W2 (111) 


2005 


Dec 


05 


13 
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Table 3. Calibrator Star Angular Diameters 



Calibrator 




^ogg 


Spectral 


E(B - V) 






Name 


(K) 


(cm s"^) 


Classification 


(mag) 


(mas) 


Ref. 


HD 6210 


6087 


3.73 


F6 V 


0.00 


0.523 ±0.027 


1, 2, 3, 4 


HD 6961 


7762 


3.80 


A7 V var 


0.00 


0.608 ±0.019 


2 


HD 32977 


8511 


4.19 


A5 V 


0.00 


0.344 ± 0.033 


2 


HD 43042 


6556 


4.28 


F6 V 


0.00 


0.590 ±0.016 


3 


HD 107193 


8710 


3.93 


AO Vn 


0.04 


0.308 ± 0.005 


2 



References. — 1. Cayrel de Strobel et al. (1997); 2. AUende Prieto & Lambert (1999); 
3. Lambert & Reddy (2004); 4. Nordstrom et al. (2004). 
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Table 4. Calibrated K' Visibilities 



Target 


Date 


Orbital 


Telescope 


Baseline 


-I— \ A 

ros. Ang. 




Name 


(HJD-2,400,000) 


Phase 


Pair 


(m) 


/I \ 

(deg) 


Visibility'' 


^ (If* CI 


cjoqi q 701 


1 97 


W1 /S9 

V V X / OZi 


947 37 


345 98 


n 445 + 016 


y wCto 


■iPQI Q 731 


1 97 


W1 /S9 

V V ± / kJZi 


947 OQ 


349 34 

OrtZ . xJrt 


n 448 + 016 


y v^ctG 


i^oqi q 749 


1 97 


W1 /S9 

V V X / OZi 


946 71 


339 99 


AQO -1-0 0-17 






1 98 


Wl /S9 

VV X / kJZ 


946 98 


336 45 


-1- 017 


0/ 1 , Q 

y v^do 


cioqi q 7c;q 


198 


Wl /S9 

VV X / kJZ 


945 87 


'^'^4 14 

OOt:. XH: 


fi08 4- 019 

U.UUO -J— U.UXc/ 


J Kjdo 


^9Q99 7Q4 


n 149 


Wl /S9 

V V X / kJZ 


949 '^l 

^rtZj. xj X 


S91 8^^ 
oz X . 00 


4qR -u 01 6 


0/ 1 , a G 

y v^cLo 


^9Q99 Sn4 


n 1 4*^ 


Wl /S9 

VV X / kJZ 


941 19 

Zjrt X . XZ 


SI Q 06 

X C/ .uu 


c;99 -u 01 6 


y v^do 


'i9Q99 81 9 


14S 


Wl /S9 

V V X / OZj 


940 06 


'^l 6 88 


KQK -1-0 015 
yj.Oc/O ZLL u.uxo 


0/ (If* Q 

y v-^do 




143 


Wl /S9 

VV X / kJZi 


93Q 06 


31 5 01 

oxu .ux 


c:qQ -1- 018 

U.Ol/O ZLl U.UXO 


^ (If* CI 

y v^do 




14"^ 


Wl /S9 

V V X / OZi 


9'^8 1 1 

ZOO. X X 


■^1 36 

0x0. ou 


c;q9 -1- 094 

U.Ot/Z 1 U.UZrt 


y wdo 


'i9Q99 834 


143 


Wl /S9 

V V X / kJZi 


936 47 


31 77 

OXU. / / 


f;94 -1- 018 

u.uzrt 1 u.uxo 


J' v^db 


'i9q9fi 747 


1 69 


Wl /<^9 

VV X / oz 




009 09 


4c;« -L 019 


0/ 1 Q C 


^9q9fi 71^4 


n 1 fi9 


Wl /S9 

VV 1 / OZi 




oou . x^ 


QS1 -u 014 




^9Q9fi 7f,"{ 


1 69 


Wl /S9 

VV X / kJZ 


944 98 

Ztit:. ZO 


"{97 ^4 

oz 1 .Ot: 


494 + 018 


'*\/ C 

J v_ d o 


^9Q9fi 77n 


1 69 


Wl /S9 

V V X / kJZ 


943 70 


395 68 

OZO . VJO 


509 + 016 

U.OU<_/ _l_ u.uxu 


1 V^do 


^9q9fi 777 


n 1 fi9 


Wl /S9 
VV 1 / oz, 


94^? 


'^9'? 78 


n p:oo I n ni 
u.ooo ZL u.uxo 


I' V^do 


^9q98 T'?! 


n 1 79 

U . -L / Z/ 


Wl /'^9 

VV X / kjZ/ 


94fi 08 


98 


n 901 _|_ n n9R 
u.zox zc u.uzu 


0/ (Ifici 
y wdo 


i^oqoQ 740 


1 79 


Wl /S9 

V V X / kJZi 


945 57 


■^■^9 66 
ooz .uu 


9fifi -1- 093 

u.zuu ZLi yj.yj^o 


^ (Ifici 
y wdo 


c;9qoQ 74fi 


1 79 


Wl /S9 

V V X / OZi 


945 90 

ZrtO . ZL/ 


■^■^1 04 

OOX .Ut: 


Q1Q -1- 016 

U.OXO ZLl U.UXU 


y wdo 


'i9Q98 7'i3 


1 79 


Wl /S9 

V V X / OZi 


944 64 


■^98 8"^ 
ozo.oo 


987 + 010 

U.ZO/ ZLl U.UXU 


^' do 


7fi1 


172 


Wl /S9 

V V X / kJZi 


944 00 

^ X X. L/L/ 


326 62 

OZ w . wz 


9q« -L 01 1 

U.Zt/O 1 U.UXi 


7 Cas 


52928.770 


0.172 


W1/S2 


243.21 


324.23 


0.303 + 0.009 


7 Cas 


53655.813 


0.743 


W1/W2 


107.91 


100.65 


0.826 + 0.019 


7 Cas 


53655.826 


0.743 


W1/W2 


107.91 


96.49 


0.790 + 0.018 


7 Cas 


53655.842 


0.743 


W1/W2 


107.70 


91.62 


0.791 + 0.021 


7 Cas 


53655.854 


0.743 


W1/W2 


107.34 


87.65 


0.816 ±0.019 


7 Cas 


53655.866 


0.743 


W1/W2 


106.84 


83.78 


0.807 ±0.019 


7 Cas 


53655.879 


0.743 


W1/W2 


106.19 


79.89 


0.756 ± 0.019 


7 Cas 


53655.891 


0.743 


W1/W2 


105.37 


75.74 


0.781 ± 0.017 


7 Cas 


53655.904 


0.743 


W1/W2 


104.45 


71.72 


0.777 ± 0.017 
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Table 4 — Continued 



Target 


Date 


Orbital 


Telescope 


Baseline 


T— \ A 

ros. Ang. 




Name 


(HJD-2,400,000) 


Phase 


Pair 


(m) 


(deg) 


Visibility^ 


J KJCXij 




743 


W1 /W2 

V V i / \ \ ^ 


103 41 


67 63 


n 770 -L f) niQ 

W.I I U _l_ W.WXtJ 


0/ (If* Q 




744 


W1 /W9 

V V X / V V Z( 


102 19 


63 12 

UO. XZi 


f) 709 _|_ f) niq 


J KJCXij 




744 


W1 /W2 

V V X / V V 


100 93 


58 68 


f) 7Q5 + n 018 

W. 1 KJ\J _l_ W.WXtJ 


r V_, CXO 


53655 961 


744 

\J . 1 J- J- 


W1/W2 

V V X / V V ^ 


99 00 


52 04 

xj W^ 


n 741 + n noi 

W • 1 ^ X 1 w • W ^ X 




536^6 652 


747 


Wl /W2 

V V X / V V Zi 


QQ 51 


15513 

xww. xo 


839 + 024 

\J.kJkjU _I_ W.WZjt: 




53656 664 


747 


Wl /W2 

V V X / V V Zi 


1 00 24 


1 50 39 

X w w . c/ 


809 + 023 




53656 678 


747 


Wl /W2 

V V X / V V Zi 


101 09 

X U X . uc/ 


145 40 

Xt:W . t:W 


830 + 029 


r ClO 


53656 6Q1 


747 


Wl /W2 

V V X / V V Z; 


101 97 

X W X . t7 1 


140 64 

XtiW . Wt: 


n 709 + n 027 

W . 1 tyZj _l_ W . WZj 1 




53656 707 


747 


Wl /W2 

V V X / V V 


103 12 


1 34 75 




y v^'cio 


53656 720 


747 


Wl /W2 

V V X / V V Z/ 


104 04 


1 30 07 

X(JW.W i 


n QQQ -L n 024 




53656 734 


747 


W1/W2 

V V X / V V ^ 


104 90 


1 25 58 

X^O .(JO 


n Q^q _|_ f) noi 

W.UOt-/ —1— w.w^x 


y v^do 


53656 747 


748 


Wl /W2 

V V X / V V Z; 


1 05 70 


121 17 

X^ X . X i 


W.(J(J(J _l_ W.WZ/X 




5'^656 769 


748 


Wl /W2 

V V X / V V Zj 


106 52 


1 1 6 04 

X X U . Wt: 


n 070 -L n noo 

U.O 1 U.UZjZ/ 




53656 775 


748 


Wl /W2 

V V X / V V Zi 


1 07 1 

XW 1 . xw 


1 1 1 71 

X X X . 1 X 


n 840 + n 020 

W.(Jt:W _I_ W.WZjW 




53656 788 


748 


Wl /W2 

V V X / V V Zi 


1 07 51 

X W 1 . tJ X 


1 07 71 

XW 1 . 1 X 


n 884 + 025 




53658 701 


757 


Wl /W2 

V V X / V V Zi 


1 03 1 1 

X WO . X X 


1 34 83 

X Ot: . (JO 


875 + 021 

W . (J 1 W _l_ W . WZj X 


y V^dO 


5^^658 714 


757 


Wl /W9 

V V X / V V Z( 


10*^ Q7 


1 '^O 47 

XOU .t: 1 


n oon 1 n non 

U.OOU _1_ U.UZjW 




53658 727 


757 


Wl /W2 

V V X / V V Z/ 


104 86 


1 25 83 


f) + n 018 

W.(J(JX _l_ W.WXtJ 


J KJCmj 


53658 741 


757 


Wl /W2 

V V X / V V Z; 


1 05 67 

X W(J . W i 


121 35 

X^ X .KJ\J 


n 070 -L n niq 

W.iJ 1 _l_ W.WXty 




53658 81 4 


758 


Wl /W2 

V V X / V V Z/ 


107 93 

X W 1 . UkJ 


97 77 

tJ i . i i 


f) 77q _|_ f) noi 

W. I I U _l_ W.WZ/X 




53658 827 


758 


W1/W2 

V V X / V V ^ 


107 82 

X w 1 . 0^ 


93 74 

<_/0 . 1 t: 


742 + 021 

W. 1 _I W.W^X 


7 Cas 


53658.840 


0.758 


W1/W2 


107.54 


89.60 


0.740 ± 0.023 


7 Cas 


53658.854 


0.758 


W1/W2 


107.01 


84.98 


0.704 ± 0.024 


7 Cas 


53658.876 


0.758 


W1/W2 


105.88 


78.23 


0.684 ±0.022 


7 Cas 


53658.892 


0.758 


W1/W2 


104.77 


73.07 


0.740 ± 0.022 


7 Cas 


53658.906 


0.758 


W1/W2 


103.61 


68.38 


0.754 ± 0.025 


Per 


52915.807 


0.794 


W1/S2 


249.35 


333.52 


0.533 ± 0.019 


Per 


52915.818 


0.794 


W1/S2 


249.21 


330.46 


0.520 ±0.016 


(f) Per 


52915.829 


0.794 


W1/S2 


248.93 


327.43 


0.475 ± 0.015 


Per 


52915.840 


0.794 


W1/S2 


248.50 


324.62 


0.540 ± 0.016 



Table 4 — Continued 



Target Date Orbital Telescope Baseline Pes. Ang. 



Name 


(HJD-2,400,000) 


Phase 


Pair 


fm 

V 


) 


(deg) 




Visibility^ 


Per 


52915.850 


0, 


.794 


W1/S2 


247, 


.94 


322. 


,06 


0, 


.586 ± 0.020 


(f) Per 


53655.657 


0, 


.635 


W1/W2 


89, 


.21 


165. 


,54 


0. 


,771 ± 0.028 


Per 


53655.673 


0, 


.635 


W1/W2 


90, 


.30 


159. 


,13 


0. 


,753 ± 0.022 


(f) Per 


53655.687 


0, 


.635 


W1/W2 


91, 


.70 


153, 


,00 


0. 


,783 ±0.027 


(f) Per 


53655.702 


0, 


.635 


W1/W2 


93, 


.31 


147. 


,18 


0. 


,743 ±0.020 


(f) Per 


53655.715 


0, 


.635 


W1/W2 


94, 


.94 


142. 


,08 


0, 


,789 ±0.021 


Per 


53655.728 


0, 


.635 


W1/W2 


96, 


.63 


137. 


,21 


0. 


,724 ±0.019 


(j) Per 


53655.743 


0, 


.635 


W1/W2 


98, 


.54 


132. 


,02 


0. 


,749 ±0.020 


Per 


53655.758 


0, 


.635 


W1/W2 


100, 


.43 


127. 


,03 


0. 


,753 ± 0.018 


(f) Per 


53655.772 


0, 


.635 


W1/W2 


102, 


.16 


122. 


,38 


0. 


,761 ±0.019 


Per 


53655.789 


0, 


.636 


W1/W2 


104, 


.03 


117. 


,05 


0. 


,724 ± 0.018 


Per 


53656.864 


0, 


.644 


W1/W2 


107, 


.89 


94. 


,57 


0. 


,836 ± 0.025 


Per 


53656.878 


0, 


.644 


W1/W2 


107, 


.55 


90. 


,58 


0. 


,866 ±0.024 


Per 


53656.892 


0, 


.644 


W1/W2 


106, 


.92 


86. 


,85 


0, 


,840 ± 0.024 


Per 


53656.905 


0, 


.644 


W1/W2 


106, 


.00 


83, 


,15 


0. 


,827 ±0.024 


Per 


53656.932 


0, 


.645 


W1/W2 


103, 


.25 


75, 


,29 


0, 


,820 ±0.031 


Per 


53656.947 


0, 


.645 


W1/W2 


101, 


.27 


70. 


,73 


0. 


,875 ±0.026 


Per 


53656.960 


0, 


.645 


W1/W2 


99, 


.35 


66. 


,67 


0. 


,887 ± 0.026 


Per 


53657.690 


0, 


.651 


W1/W2 


92, 


.54 


149. 


,83 


0. 


,937 ± 0.039 


Per 


53657.703 


0, 


.651 


W1/W2 


94, 


.10 


144. 


,64 


0. 


,880 ± 0.031 


Per 


53657.717 


0, 


.651 


W1/W2 


95, 


.89 


139, 


,29 


0. 


,744 ±0.024 


Per 


53657.733 


0, 


.651 


W1/W2 


97, 


.90 


133. 


,73 


0. 


,684 ±0.022 


Per 


53657.746 


0, 


.651 


W1/W2 


99, 


.59 


129. 


,25 


0. 


,643 ±0.021 


Per 


53657.759 


0, 


.651 


W1/W2 


101, 


.25 


124. 


,84 


0. 


,554 ±0.018 


Per 


53657.775 


0, 


.651 


W1/W2 


103, 


.18 


119. 


,53 


0. 


,620 ±0.018 


Per 


53657.790 


0, 


.651 


W1/W2 


104, 


.65 


115. 


,13 


0. 


,707 ±0.021 


Per 


53657.879 


0, 


.652 


W1/W2 


107, 


.40 


89. 


,50 


0. 


,725 ± 0.027 


Per 


53657.892 


0, 


.652 


W1/W2 


106, 


.71 


85. 


,92 


0. 


,826 ± 0.020 


Per 


53657.905 


0, 


.652 


W1/W2 


105, 


.77 


82. 


,34 


0. 


,861 ± 0.021 


Per 


53657.925 


0, 


.652 


W1/W2 


103, 


.73 


76. 


,48 


0. 


,924 ± 0.017 
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Table 4 — Continued 



Target 


Date 


Orbital 


Telescope 


Baseline 


T— \ A 

ros. Ang. 




Name 


(HJD-2,400,000) 


Phase 


Pair 


(m) 


(deg) 


Visibility^ 


\y ± Kji. 




653 


W1 /W2 

V V i / y\ £j 


102 16 


72 72 

1 1 Zj 


n oQi + n 016 

W.(-/OX _l_ w.wxw 


/Tl PpT* 


UOUU / .C/UVJ 


653 


W1 /W2 

V V X / V V Z( 


100 36 


68 77 


f) Qoq -1- n 016 


t Tan 




242 


SI /R1 


329 84 


33 00 

OO. WW 


f) 004 + n 037 

yj.u^rt _l_ W.WOi 


C Tan 


53360 774 


0.242 


SI /El 


328 55 


30 86 

o w . ow 


W.iyXO —1— W.WO^ 


C Tan 


5336(1 7Q5 


242 


SI /El 

kJ X / XJ X 


325 81 


27 32 


898 + 036 

W.(Jty(J _l_ w.wow 


Tan 


533fin 807 


242 


SI /El 

kJ X / XJ X 


323 78 


24 94 


789 + 038 

W. 1 OC/ _l_ w.woo 


Tail 


533fin 845 


242 


SI /El 

VJ j XJ X 


31 7 39 

OX I .xJO 


1 7 02 

X 1 . WZi 


615 + 021 

W.WXO —1— W.WZiX 


Tan 

^ -L CI L L 


533fin 850 


242 


SI /El 


31 6 65 


1 5 94 

X O . i^t: 


n Ran + n noo 

W.WWty _l_ W.WZjW 


Tan 


533fi1 7fin 


24Q 


SI /El 

OX / XjX 


329 62 


32 58 


W.OXO ZLl U.UOO 


C Tan 


534fi1 fifi2 


001 


W1 /SI 

V V X / kJ X 


261 82 


304 56 

OWrt .o w 


W.OWO _l_ W.WXty 


Tan 


534fi1 fi73 


001 

\J .\J\J ±. 


W1 /SI 

V V X / kJ X 


256 24 


303 94 

OWO. Urt 


f) 400 _|_ f) niq 

W.rtOL? _l_ W.WXty 


C Tan 


53708 7Q5 


860 


W1 /W2 

V V X / V V 


92 98 


108 47 


n 74Q _|_ n 024 

W. i rtU _l_ W.W.^t: 


t Tan 


K,Q7n8 8nQ 

OO 1 UO. OUiy 


860 


W1 /W2 

V V X / V V Z( 


96 78 


106 04 

XWU. Wt: 


n yiq _|_ n noQ 

W. 1 Xc/ W.WZjO 


Tan 


53708 822 


860 


Wl /W2 

V V X / V V Zi 


1 00 05 


103 84 

X WO . (Jt: 


803 + 022 

W.(JWO _l_ W.WZjZj 


Tan 


53708 838 


861 


Wl /W2 

V V X / V V Zi 


1 03 33 

X WO . oo 


101 35 

X W X . OO 


775 + 028 

W. 1 1 _l_ W.WZiO 


C Tan 

L J. fX Li 


53708 853 


861 


Wl /W2 

V V X / V V Zi 


105 68 

X WO . wo 


99 12 

U U . XZj 


746 + 022 

W. ( ^W —1— W.WZjZj 


t Tan 


5^^708 8fiQ 

1 UO. OUc/ 


861 
u .ou ± 


Wl /W2 

V V X / V V Z( 


107 24 

X W 1 . 


96 97 


W. 1 U 1 _I_ W.UZ/Z/ 


C Tan 


53708 87fi 


861 


Wl /W2 

V V X / V V 


107 67 

X W 1 .\J I 


95 97 


n 7-]« -1- n noi 

W. 1 XW _l_ w.w.^x 


Tan 


53708 8Q0 

XJKJ I \J\J .yDUKJ 


861 


Wl /W2 

V V X / V V 


107 93 

X W 1 . Ukj 


94 14 


n a-jFt + 020 

W. W 1 _l_ w. w.^w 


C Tan 


5370Q 71 8 


867 


Wl /W2 

V V X / V V 


65 27 


129 49 


n 740 _i_ n nqo 


C Tan 


53700 734 


867 


W1/W2 

V V X / V V Zj 


71 89 


123.24 


924 + 029 

W . _I W • W Zj <_/ 


CTau 


53709.751 


0.867 


W1/W2 


78.20 


118.27 


0.868 ±0.028 


C Tau 


53709.765 


0.868 


W1/W2 


83.77 


114.39 


0.946 ±0.029 


C Tau 


53709.782 


0.868 


W1/W2 


89.40 


110.74 


0.927 ±0.028 


C Tau 


53709.795 


0.868 


W1/W2 


93.81 


107.95 


0.948 ± 0.032 


CTau 


53709.809 


0.868 


W1/W2 


97.51 


105.57 


0.815 ±0.021 


CTau 


53709.823 


0.868 


W1/W2 


100.87 


103.25 


0.803 ± 0.029 


CTau 


53709.840 


0.868 


W1/W2 


104.15 


100.64 


0.720 ± 0.019 


K Dra 


53462.655 


0.074 


S1/S2 


26.40 


29.89 


0.929 ± 0.037 


K Dra 


53462.756 


0.076 


S1/S2 


27.74 


4.61 


0.995 ± 0.021 
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Table 4 — Continued 



Target 


Date 


Orbital 


Telescope 


Baseline 


ros. Ang. 




Name 


(HJD-2,400,000) 


Phase 


Pair 


(m) 


(deg) 


Visibility^ 


ixj J_/1 CX 


53462 774 


076 


SI /S2 


27 77 


27 




rv J_/i d 




077 


SI /S2 


27 70 


■^53 51 


1 nno _|_ f) noK 




53462 827 


077 

\j .yj 1 1 


SI /S2 


27 51 


346 82 


1 QKQ -L n 041 


rv J— /± d 


53462 845 


077 


S1/S2 


27 29 


342.31 


990 + 023 




53462 S5Q 


078 


SI /S2 

kJ -L / kJ.^ 


27 09 


338 86 

^J^J(_y . (J VJ 


1 01 9 + 034 




53462 QDQ 


078 


SI /S2 


26 04 


326 47 


1 018 + 028 

X.WXO —1— w.w^o 




53462 Q26 


079 


SI /S2 


25 55 


322 1 3 


974 + 035 


rv J— '± Oj 


53462 Q46 


079 


SI /S2 

kJ X / kJ.^ 


24 92 


31 7 1 8 

OX 1 . X o 


1 (-127 + n 025 


LT TlTO 


53462 Q63 


079 


SI /S2 

kJ± / OZi 


24 27 


31 2 66 


070 _|_ n noi 


/<• T)t'Q 

My J_/l d 


53462 Q82 


080 


SI /S2 

kJ i / kJ.tj 


23 50 


307 78 


1 OAK -L n 024 




53462 QQQ 


080 


S1/S2 

kJ X / kJ^ 


22 79 


303 56 


n nqq 4- n qao 

yj .tj tjtj —1— yj .yj^^ 


My J_/± d 


53466 750 


0.141 


Sl/El 

kJ X / J-Ji 


267 57 


43 26 


n «Q« _|_ n 022 


My J_/i d 


5'^4fifi 847 


142 


SI /F,1 

kJ X / 11/X 


28(S 48 


1 95 

XU . C/tJ 


n 7QQ -L n nQ2 


AC TjT'i^ 


53466 87(1 


143 


SI /R1 

kJ X / XJ X 


288 46 


1 67 


727 + 032 


tC 1 IT 


5S4fifi 8Qn 


14'^ 


SI /F,1 

kJ X / J_/X 


289 41 


5 '^7 


n 714 _)_ n QQi 

U. 1 X^ _1_ u.uox 


AC TjTfl 
rb J^i d 


53466 Ql 4 


144 


SI /R1 

kJ X / XJ X 


289 72 


358 90 


681 + 033 

W.VJOX _1_ \J .\JxJkJ 


i- - T ITQ 

rb J— /i d 


c,Q47'? 794 


254 


SI /F,1 

kJ ± / J_J J- 


2(i5 84 


AA 8(S 

tit: . OU 


n 77« _L n noq 


AC TlTP* 

My J_/l d 


53708 QQ5 


076 

yj .yj 1 yj 


W1 /W2 

V V i / V V 


104 62 


1 54 45 


f) CQ2 + n 027 

yj.KJKj^ _i_ yj.yj^i 


AC T)!*?* 
My J_/l d 


537nQ (1(18 


077 

yj .yj 1 1 


W1 /W2 

V V i / V V 


1 04 95 


149 84 


n Q 55 -L n 057 


AC TjTi^ 
My J_/l d 


537(1Q 020 


077 

yj .yj 1 1 


W1 /W2 

V V i / V V 


1 05 29 


1 45 59 


f) Q47 _|_ f) 


My J— /± d 


5370Q 868 


090 


W1/W2 

V V X / V V ^ 


104.32 


200 84 


817 + 026 


K Dra 


53709.883 


0.091 


W1/W2 


104.05 


195.44 


0.713 ±0.018 


K Dra 


53709.897 


0.091 


W1/W2 


103.85 


190.28 


0.753 ±0.023 


K Dra 


53709.910 


0.091 


W1/W2 


103.73 


185.24 


0.859 ±0.021 


n Dra 


53709.925 


0.091 


W1/W2 


103.69 


179.64 


0.835 ±0.018 


K Dra 


53709.937 


0.092 


W1/W2 


103.73 


175.03 


0.760 ± 0.017 


K Dra 


53709.950 


0.092 


W1/W2 


103.83 


170.33 


0.719 ±0.015 


K Dra 


53709.962 


0.092 


W1/W2 


104.00 


165.71 


0.799 ± 0.018 


K Dra 


53709.975 


0.092 


W1/W2 


104.21 


161.15 


0.873 ± 0.015 


K Dra 


53709.987 


0.092 


W1/W2 


104.48 


156.65 


0.889 ± 0.018 
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Table 4 — Continued 



Target 
Name 


Date 
(HJD-2,400,000) 


Orbital 
Phase 


Telescope 
Pair 


Baseline 

(m) 


Pos. Ang. 

(deg) 


Visibility^ 


K Dra 

K, Dra 
K Dra 


53710.000 
53710.031 
53710.043 


0.093 
0.093 
0.093 


W1/W2 
W1/W2 
W1/W2 


104.80 
105.70 
106.07 


151.94 
140.57 
136.19 


0.926 ±0.022 
0.879 ±0.026 
0.871 ± 0.029 



'^The formal visibility errors listed here may underestimate the actual errors by a factor of 
«i 2.8 (§2.2). 



Table 5. Spectroscopic Observations 







Resolving 


Range 


Date 




Telescope 


Spectrograph 


Power 


(A) 


(UT) 


Targets 


KPNO coude feed 0.9 m 


Coude 


12800 


4235-4586 


2004 Aug 16-21 


7 Cas, Per 


KPNO coudc feed 0.9 m 


Coude 


12800 


4236-4587 


2005 Nov 10-14 


7 Cas, (f) Per, ( Tau, k Dra 


KPNO coude feed 0.9 m 


Coude 


9150 


6465-7140 


2004 Oct 12-14 


7 Cas, (j) Per, ( Tau 


McDonald Obs. 2.7 m 


CoolSpec 


1520 


20691-22520 


2004 Feb 06 


7 Cas, 4> Per, ( Tau 


McDonald Obs. 2.7 m 


LCS 


810 


3020-5717 


2005 Apr 29 


K Dra 


Ondfejov Obs. 2.0 m 


Coude 


12800 


6262-6774 


2005 Apr 11 


K Dra 


Ritter Obs. 1.0 m 


Echelle 


26000 


5413-6596 


2005 Apr 10 


K Dra 



Table 6. Thick Disk Model Fits of Visibilities 





7 Cas 


7 Cas 


(p Per 


(p Per 


C Tau 


C Tau 


K Dra 


K Dra 


Parameter 


Single 


Binary 


Single 


Binary 


Single 


Binary 


Single 


Binary 


i (deg) 


51 ±4 


127 ±4 


90 ±7 


69 ±5 


90 ±3 


89 ±2 


72± 18 


26 ±9 


a (deg) 


116 ±4 


114 ±3 


46 ±4 


49 ±3 


38 ±2 


37 ±2 


114 ±4 


21 ±3 


logpo fit (g cm-3) ... 


-10.14 ±0.24 


-10.18 ±0.25 


-10.92 ± 0.10 


-10.98 ±0.08 


-9.71 ± 0.31 


-9.73 ± 0.35 


-12.55 ± 0.39 


-12.21 ±0.32 


logpo IRAS* (g cm-3) 


-10.95 ± 0.35 


-10.95 ±0.35 


-10.85 ± 0.25 


-10.85 ±0.25 


-11.15 ± 0.25 


-11.15 ±0.25 


-11.35 ± 0.15 


-11.35 ±0.15 


"fit 


2.70 ±0.31 


2.65 ±0.32 


1.80 ± 0.09 


1.76 ±0.08 


3.14 ±0.38 


3.19 ±0.40 


0.16 ±0.43 


0.67 ±0.36 


niRAS* 


3.25 


3.25 


3.0 


3.0 


3.25 


3.25 


3.0 


3.0 




24.3 


24.0 


8.0 


7.1 


7.9 


7.9 


4.0 


4.0 


E(V - K)fit (mag) . . . 


1.60 


1.53 


0.58 


0.63 


0.53 


0.40 


0.31 


0.14 


E(V - KU,'^ (mag) . 


0.85 


0.85 


0.68 


0.68 


0.65 


0.65 


0.39 


0.39 


Old (mas) 


0.50 


0.50 


0.30 


0.30 


0.40 


0.40 


0.39 


0.39 


Disk FWHM (mas) . . 


1.35 


1.36 


0.98 


0.51 


1.99 


1.82 


3.02 


1.72 



^'Waters et al. (1987). 
''Dougherty et al. (1994). 
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Table 7. Gaussian Elliptical Fits of Visibilities 



Parameter 


7 Gas 


(f) Per 


CTau 


K Dra 


6'gd fit (mas) . 


1.95 ±0.07 


2.30 ±0.08 


1.79 ±0.07 


1.83 ±0.11 


Ogd Ua"^ (mas) 


3.47 ±0.02 


2.67 ±0.20 


4.53 ±0.52 




Ogd Ha^ (mas) 


3.59 ±0.04 


2.89 ±0.09 


3.14 ±0.21 




^^fit 


0.59 ± 0.04 


0.00 ±0.22 


0.09 ±0.22 


0.00 ± 0.67 


rna^ 


0.70 ± 0.02 


0.46 ± 0.04 


0.28 ±0.02 




THa'' 


0.58 ± 0.03 


0.27 ±0.01 


0.31 ± 0.07 




0fit (deg) .... 


25 ±4 


-44 ±3 


-52.2 ±1.7 


-8 ±9 


0Ha''(dcg) ... 


19 ±2 


-62 ±5 


-58 ±4 




0Ha''(deg) ... 


31.2± 1.2 


-61.5 ±0.6 


-62.3 ±4.4 




Cp fit 


0.395 ±0.018 


0.538 ±0.015 


0.414 ±0.029 


0.628 ± 0.030 




23.6 


7.9 


7.7 


8.6 



^Quirrenbach et al. (1997). 
''lycner et al. (2004, 2006). 



-41 - 



Table 8. Emission Line Equivalent Widths and Density Fits 



Parameter 


7 Cas 


(p Per 


( Tau 


K Dra 


W^A(Ha) (A) 


-31.2 


-42.6 


-19.5 


-20.3 


1^a(H7) (A) 


+1.0 


+2.1 


+4.0 


+5.9 


W^A(Br7) (A) 


-5.8 


-9.2 


-2.8 




logpo fit (Ha) (g cm-3) 


-11.7 


-12.5 


-11.2 


-13.4 


logpo fit(H7) (g cm-3) 


-11.7 


-12.8 


-11.2 


-13.7 


logpo fit(Br7) (g cm~^) 


-12.0 


-12.9 


-11.7 




logpo &t{K') (g cm-3) . 


-10.1 ±0.2 


-11.0 + 0.1 


-9.7 + 0.3 


-12.2 + 0.3 
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Fig. 1. — The left panel is a depiction of the K' image of a disk model for the case of 7 Cas 
[i = 51°, po = 7.2 X 10~^^ g cm^^, n = 2.7). The right panel shows the corresponding 
interferometric visibility V as a function of projected baseline for position angles along the 
apparent minor axis {top line) and major axis {bottom line). 



1.0 0.5 0.0 -0.5 
(mas) 



-1.0 



50 100 150 200 250 300 350 



Fig. 2. — An image and visibility plots for a 7 Cas model as in Fig. 1 but with a larger 
inclination oi i — 80°. The dotted hnes in the visibihty panel repeat the curves from the 
original model in Fig. 1. 
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Fig. 3. — An image and visibility plots for a 7 Cas model as in Fig. 1 but with a smaller 
base density of po — 3.6 x 10~^^ g cm~^. The dotted hnes in the visibility panel repeat the 
curves from the original model in Fig. 1. 
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Fig. 4. — An image and visibility plots for a 7 Cas model as in Fig. 1 but with a smaller 
radial density exponent of n = 2.0. The dotted lines in the visibility panel repeat the curves 
from the original model in Fig. 1. 




Fig. 5. — The observed Ha line in 7 Cas {thick line) compared with the synthetic profile 
for a Keplerian disk {double-peaked thin line) and the same convolved with a Voigt profile 
{broader thin line). 
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Fig. 6. — The disk model for 7 Cas. The left panel shows the appearance of the model star 
and disk in the sky (with the spin axis rotated east from north by angle a) . The right panel 
shows the associated Fourier transform of the spatial image, which is directly related to the 
predicted interferometric visibility. The plus signs in the upper half of the panel are scaled 
representations of the observed visibility squared at the {u, v) coordinate of observation, 
while the mirror symmetric, thick line segments in the lower half show the scaled residuals 
from the fit (vertical for a positive residual and horizontal for a negative residual). 
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Fig. 7. — Hydrogen line profiles in the spectrum of 7 Cas plotted as a function of heliocentric 
radial velocity. The solid, thick lines indicate spectra made in 2004 while the dashed, thick 
lines correspond to 2005 observations. The thin lines represent the model synthetic spectra. 
The continuum levels for H« (top) and H7 {middle) are offset in rectified flux by +2 and 
+ 1, respectively, for clarity of presentation. 




Fig. 8. — The model image and visibilities for the binary star fit of the interferometric data 
of (f) Per in the same format as Fig. 6. The secondary appears as the bright dot in the lower 
part of the spatial image. 




Fig. 9. — Hydrogen line profiles in the spectrum of (j) Per plotted in the same format as 
Fig. 7. 
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Fig. 10. — The model image and visibilities for the single star fit of the interferometric data 
of ( Tau in the same format as Fig. 6. 




Fig. 11. — Hydrogen line profiles in the spectrum of ( Tau plotted in the same format as 
Fig. 7, except that the Br7 profile is offset by —0.5 in flux for clarity. 
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Fig. 12. — The model image and visibilities for the binary star fit of the interferometric data 
of K Dra in the same format as Fig. 6. The secondary appears as the bright dot in the upper 
left part of the spatial image. 
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Fig. 13. — Hydrogen line profiles in the spectrum of k Dra plotted in the same format as 
Fig. 7. We lack observations of Br7 for this star but we include our low resolution spectrum 
of H/3 instead (without a synthetic comparison). The continuum levels for H« {top) and B.f3 
{middle) are offset in rectified flux by +2 and +1, respectively, for clarity of presentation. 
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Fig. 14. — The summed spatial intensities projected onto the major axis for our single star 
model of 7 Cas. The solid line shows the K'-hand summed intensity while the dashed line 
shows the summed intensity for an Ha image integrated over a 2.8 nm bandwidth. Both are 
renormalized for convenience to their respective summed intensity at i? = 1.14i?^. 



